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ABSTRACT 

o : 
o . 

CN ■ The near-IR and X-ray flares in Sagittarius A* are believed to be produced 

Ph 1 by relativistic electrons via synchrotron and synchrotron self-Comptonization, 

respectively. These electrons are likely energized by turbulent plasma waves 

\Q • through second order Fermi acceleration that, in combination with the radiative 

cooling processes, produces a relativistic Maxwellian distribution in the steady 

CN ■ state. This model has four principal parameters, namely the magnetic field B, 

\q • the electron density n and temperature 7 c m e c 2 , and the size of the flare region 

R. In the context of stochastic acceleration, the quantities Rn x ' 2 B and ^ c Rn 

fC) ■ should remain nearly constant in time. Therefore, simultaneous spectroscopic 

observations in the NIR and X-ray bands can readily test the model, which, 
if proven to be valid, may be used to determine the evolution of the plasma 

r£3 ■ properties during an eruptive event with spectroscopic observations in either band 

or simultaneous flux density measurements in both bands. The formulae can be 
applied to other isolated or confined systems, where electrons are accelerated to 
relativistic energies by plasma wave turbulence and produce most of the emission 
!• ■ via synchrotron processes. 
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1. Introduction 

Sagittarius A*, the compact radio source at the Galactic Center, is powered by accretion 
onto an M ~ 3-4 x 10 6 M© supermassive black hole (Schodel et al. 2002; Ghez et al. 2004; 
Melia 2006). There is compelling evidence that the near-IR (NIR) flares in Sagittarius A* are 
produced by relativistic electrons via synchrotron emission within a few Schwarzschild radii 
of the black hole (Genzel et al. 2003; the Schwarzschild radius r s = 10 12 (M/3.4 x 10 6 M ) 
cm) . Inverse Compton scattering of the NIR photons by the same electrons produces X- and 
7-ray emission, which may account for the observed X-ray flares (Baganoff et al. 2001, 2003; 
Goldwurm et al. 2003; Porquet et al. 2003; Belanger et al. 2006). Simultaneous multi- 
wavelength observations during a handful of flares have generally supported this hypothesis 
(Eckart et al. 2004, 2005; Belanger et al. 2005; Yusef-Zadeh et al. 2006). However, 
our knowledge of the flare energization and electron acceleration mechanisms is still very 
rudimentary, though some details are starting to emerge (Liu et al. 2004, 2006; Tagger & 
Melia 2006; Bromley, Melia, & Liu 2001; Broderick & Loeb 2005). 

The first NIR spectroscopic observations at a frequency v ~ 1.4 x 10 14 Hz revealed very 
soft power spectra, vF v oc v a , with index —4 < a < — 2 (Eisenhauer et al. 2005), indicating a 
sharp cutoff in the electron distribution. Such a cutoff is a natural consequence of stochastic 
acceleration (SA) by plasma waves, in which low-energy electrons are accelerated up to 
a relativistic energy, where radiative cooling becomes important (Schlickeiser 1984; Park & 
Petrosian 1995). In Sagittarius A*, the X-ray flare luminosity is always lower than that of the 
NIR events, so synchrotron cooling dominates. More recent NIR observations also indicate 
that F u may be correlated with a (Ghez et al. 2005; Gillessen et al. 2006; Hornstein et al. 
2006). 

Flares in Sagittarius A* are likely driven by an accretion instability (e.g., Tagger & 
Melia 2006), which triggers the transfer of gravitational energy of the plasma into tur- 
bulence. In earlier work, we demonstrated that if thermal synchrotron and synchrotron 
self-Comptonization are the dominant emission mechanisms during Sagittarius A*'s flares, 
simultaneous X-ray and NIR spectroscopic observations (SSO) can measure the size of the 
flare region R, the electron density n and "temperature" r y c m e c 2 , and the magnetic field B 
(Liu, Melia, & Petrosian 2006). Interestingly, most of the flares are found to be local events 
with the source size R < r s , which can also explain the quasi-period modulations observed 
in some flares (see also Liu, Petrosian, & Melia 2004). This justifies the neglect of the source 
structure as a zero order approximation and suggests that the flares involve highly dynamical 
activity and the electron acceleration process has to be addressed to uncover the underlying 
physical processes. 

In this paper, we consider steady state electron populations produced via SA by tur- 
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bulent plasma waves which, in combination with the synchrotron cooling, results in a rel- 
ativists Maxwellian distribution. The electron temperature is set by the energy diffusion 
and synchrotron cooling rates. The latter depends only on B, while the former is given by 
a combination of B, n, and R, which together imply that the product 7 c ni? should not vary 
significantly as the flare evolves, nor from one flare to the next. Moreover, since electrons 
gain energy from the plasma waves, overall energy conservation requires that the luminosity 
be equal to the energy transfer rate from turbulence, which also depends only on B, n, and 
R. This requires Bn l l 2 R to be constant. These two system constraints effectively reduce 
the model parameters to just two, so that SSO can readily test the model. 

Previous work on the flares is summarized in § 2, where we discuss difference between 
these models and assess the merit of previous proposed scenarios on the definiteness of their 
predictions and/or whether the model parameters can be well constrained by observations. 
The electron acceleration is described in § 3. And in § 4 the model is applied to an NIR flare 
with a very hard spectrum. The model predictions and its limitations are discussed in § 5. 

2. Previous Work 

Several models have been proposed since the Chandra discovery of a bright X-ray flare 
from Sagittarius A* (Baganoff et al. 2001; Markoff et al. 2001; Liu & Melia 2002; Yuan et 
al. 2004; Nayakshin, Cuadra, & Sunyaev 2004). The facts that the flare variation time scale 
is comparable to the orbital period at the last stable orbit and there is a good correlation 
between the X-ray and the strongly polarized NIR emission suggests that the flares are 
produced very close to the black hole. The nature of the correlation between an X-ray flare 
and a delayed radio outburst, on the other hand, indicates that the radio emission may 
be produced at large radii (Zhao et al. 2004). In the jet model studied by Markoff et al. 
(2001), it was assumed that the electrons and protons are Maxwellian and have the same 
temperature, the plasma near the base of the jet, which is within a few Schwarzschild radii 
of the black hole, are heated suddenly to a higher temperature during flares via shocks or 
magnetic field reconnection. By changing the physical conditions at the base of the jet, 
they were able to reproduce the observed X-ray flux. The distinct spectral and polarization 
characteristics of sub-millimeter emission, on the other hand, suggests that this radiation is 
produced by a small hot magnetized accretion torus near the black hole (Melia et al. 2000, 
2001). The flat spectrum of the X-ray flares then motivated us to consider the production of 
X-rays via bremsstrahlung during an accretion instability (Liu & Melia 2002). To reduce the 
model parameters, the electrons and protons were assumed in thermal equilibrium in this 
calculation. Yuan et al. (2004) later considered a more complicated accretion model, where 
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the electrons and protons have different temperatures and there is a non-thermal component 
in the electron distribution prescribed by the choice of particle acceleration processes. With 
the temperature for the thermal component, the spectral index, high energy cutoff, and 
normalization factor characterizing the non-thermal component as free parameters, it can 
accommodate a variety of flare spectra. 

All these models are still consistent with the limited NIR and X-ray observations, due 
in part to the large number of model parameters and the relative paucity of data. It is clear 
that one must address the electron acceleration mechanism in order to (possibly) distinguish 
between these various models. Given the physical conditions in the emission region — those 
of a hot, magnetized plasma — it is natural to suppose that the electrons are being accel- 
erated by plasma waves (e.g., Dermer, Miller, & Li 1996), whose turbulent energy is fed 
from an accretion-induced instability. However, given the uncertainties in the generation of 
turbulence, and in the cascade and damping processes, several parameters were introduced 
to describe the turbulence spectrum (Liu et al. 2004). Not surprisingly, as was the case 
for the previous models, the available observations could not provide rigorous constraints on 
the parameters. Nevertheless, this first attempt at building a physical model suggested that 
some of the flares are actually produced within a region smaller than the size of the black 
hole, indicating that they are likely events localized in certain portions of the disk. 

The discovery of NIR flares with very soft spectra implies a sharp cutoff in the relativistic 
electron distribution producing the polarized NIR emission via synchrotron. For typical 
conditions near the black hole inferred from the modeling of Sagittarius A*'s quiescent 
emission (Melia et al. 2001; Yuan et al. 2003), the synchrotron cooling time of electrons 
producing NIR emission is also comparable to the flare variation time scale and may therefore 
be the cause of the cutoff. Subsequent work with this model showed that the electron 
acceleration may be insensitive to the details of the turbulence spectrum and that such a 
cutoff is a natural consequence of SA of electrons by turbulence (Liu et al. 2006). Moreover, 
if the injection of particles into the acceleration region is slow, the electron spectrum in 
steady state may be approximated as a relativistic Maxwellian distribution. In this case, 
there are only four model parameters, namely the source size R, the electron temperature 
7 c m e c 2 and density n, and the magnetic field B. This model has therefore now matured to 
the point where its parameters may be readily determined by SSO. Of course, this does not 
yet mean that SA is solely responsible for accelerating the electrons. Other mechanisms may 
also produce similar distributions. 
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3. SA of Electrons by Plasma Waves 

Tagger and Melia (2006) first showed that the Rossby wave instability in a small accre- 
tion torus can reproduce the characteristics of the flare light curves. In such a picture the 
gravitational energy of protons and ions is released by this macroscopic instability, which 
generates the turbulence on a scale comparable to R. The turbulence then cascades toward 
small scales and accelerates electrons in the process. It is clear that the whole plasma in 
the turbulent region can be energized by the gravitational energy release. We therefore do 
not need a continuous injection of electrons at low energies. Moreover, since the plasma is 
gravitationally trapped by the black hole, the accelerated electrons will not escape from the 
acceleration region as long as its energy is lower than the gravitational banding energy of the 
protons, which is comparable to the protons rest mass energy. The evolution of an electron 
distribution (7, t) under the influence of a turbulent magnetic field is given therefore by 
(Blandford & Eichler 1987) 



dN _ d 
dt c?7 



7 4 dT 2 N , 



—N 



(97 T 



(1) 



where 7 = 1007 2 is the Lorentz factor 1 and the synchrotron cooling and acceleration times 
are given, respectively, by, 



syn 



( 7 ) = r / 7 = 9m 3 e c 5 /Ae 4 B 2 -f = 21.56 B' 2 ^ 1 hrs, (2) 
r ac = 2 7 2 / < A 7 A 7 /At >= d3Rc/v 2 A = 52.49 CiR 12 n 7 B^ 2 hrs , (3) 



where v\ = B / (47mm p ) 1 ^ 2 = 6.901 x 10 8 Bin 7 ^ 2 cm s _1 is the Alfven velocity, B\ = B/10 
Gauss, n-i = n/10 7 cm -3 , i? 12 = -R/10 12 cm, and C\ is a dimensionless quantity that depends 
on the microscopic details of the wave-particle interaction, and is of order 1 if the scattering 
mean-free-path is comparable to R, which may be smaller than the Schwarzschild radius 
(~ 10 12 cm) of the black hole. The other constants have their usual meaning. Note that in 
equation (3) the velocity of the scatterers in the original SA proposed by Fermi (1949) has 
been replaced by the Alfven velocity va- This may not be always true especially when the 
sound velocity v s > va- The exact dependence of r ac on v s and va is still under investigation. 
Because v s also depends on the proton temperature, which is an input in the model, here 
we focus on the simplest scenario, where the effects of sound waves are ignored. And as 
we shall see below, particle-particle collisions can be ignored, and the electrons are mainly 



1 Here we have assumed that the direct acceleration rate < A'y/ At > /j is equal to the energy diffusion 
rate < A7A7/ At >' /j as is usually true in SA by plasma waves, where a prime indicates a derivative with 
respect to 7. In a more general case, where the direct acceleration rate < A^y/At > /-f — p < A-fAj/At >' /-f 
the first term on the right hand side of equation (1) should be [7 4p T ac 1 (7 2_4p iV)']'. 
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scattered by the turbulence. Since the inferred particle pressure is much higher than the 
magnetic field pressure, R is the only length scale in the system and should be comparable 
to the particle mean-free-path. If the flares are triggered by instabilities that are similar 
from event to event, C\ should be nearly constant. In steady state this yields 2 

iV( 7 ) = (n 7 2 /2 7c 3 ) exp(-7/ 7c ) , with 7c = r /r ac = 41.08 C^R^rij 1 . (4) 

The normalization of N is accurate to order 7 ~ 2 . Therefore 7c i?n only depends on the 
instabilities and should not change significantly with time. 



4. Relativistic Maxwellian Synchrotron Emission and Self-Comptonization 

The thermal synchrotron flux density and emission coefficient (i.e. emissivity divided 
by the solid angle) are given, respectively, by (Petrosian 1981; Mahadevan, Narayan, & Yi 
1996) 

F v = (47cR 3 /3D 2 )£ u , E v = (V3e 3 /8irm e c 2 )Bnx M I(x M ) , (5) 

where 

I( XM ) = 4.0505a; M 1/6 (l + 0.40a; M 1/4 + 0.5316a; M 1/2 )exp(-1.8899a;^ 3 ), (6) 
x M = vjv c = 47rm e cz//3e J B 7 2 = 1412 C\ is u R 2 12 n 2 B' 1 , (7) 

v = z/i4l0 14 Hz, and the distance to the Galactic Center D = D 8 8 kpc. The spectral index 
in a given narrow frequency range is therefore 

a. 1^)= 1.833- 0.6300^3 _ 01000^+0265^ 
d In v x ^ 2 + 0.4000a;^ 4 + 0.5316 

and the flux density may be written 

e 3 BnR 3 

x m I(xm) = 639.7 R\ 2 n 7 B 1 D% 2 x m I(xm) mJy. (9) 



2 v / 3m e c 2 D 2 



The thick lines in Figure 1 (left) exhibit this spectral index (dotted), and the normalized 
spectrum x 2 m I(xm) (solid), as functions of the normalized frequency 

We can estimate the energy transfer rate in a fully developed turbulence from dimen- 
sional analysis: 

£ Tur b = C 2 v A B 2 R 2 = 6.901 x 10 U C 2 R 2 12 n~ 1,2 B\ ergs s" 1 , (10) 



2 For the more general case mentioned above N(j) oc j 4p 2 exp(— 7/70). 
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where C 2 is another dimensionless quantity that only depends on the turbulence generation 
mechanism and therefore remains nearly a constant 3 . If there are no other significant energy 
loss processes, this energy transfer rate should be equal to the source luminosity (C — £ syn ) 
in the steady state, where 

£ syn = (647re 4 /9m2 c 3 ) n R 3 B 2 ^ = 8.949 x 10 U C^ 2 R 12 n^Bf ergs s" 1 . (11) 

We therefore infer that 

R12 n\ /2 B 1 = 1.297 Cr 2 ^" 1 , (12) 

a condition that should not change significantly with time either. We note that the electron 
advection occurs on the accretion time, which can be more than 10 times longer than the 
dynamical time. This term is negligible during flares (see § 5). And even though there are 
other energy loss processes, such as the diffusive escape of accelerated particles and/or the 
propagation of waves away from the flare region (Petrosian & Liu 2004), we expect that an 
approximately fixed fraction of the released energy goes into radiation. This constraint is 
still valid. 

For given values of C\ and C 2 , spectroscopic observations in the NIR band then can 
be used to determine the properties of the flaring plasma. First from equation (8) and the 
observed spectral index a, one can obtain %, which in combine with equations (7) and 
(6) and the observed frequency gives v c and I{xm)- Then one can use the measured flux 
density and equations (7), (9), and (12) to obtain R, n, and B, and 7c can be inferred 
with equation (4). For the bright NIR flare observed by Ghez et al. (2005), a = 0.5, 
z/14 = 1.429, and F u = 7 mJy, we have % = 7.314, xmI( x m) = 0.7811, v c = 1.954 x 
10 13 (xm/7.314)- 1 (z/i 4 /1.429) Hz, 

p nn o 47 i -2^6/7 n 10/7 / X M \~V 7 ( x M I(x M ) y 5/7 ( »14 \ V 7 / Fy \ ^ , , 

R 12 - 0.08471 ^ D 8 (^—j [-^7-) (jjm) [j7^) >( 13 ) 

n 7 = 2.266 C-^C-^D^ (J^f (EmIM)" (J^Y* (^-V^) 
7 128 V 7.314/ V °- 7811 / v 1-429/ V 7mJ y/ 

Bl = 10.17 C^C^D?' 7 f^L_Y 1/7 (-mI{xm)\ 217 rj^fi ( _F^_V 2/7 
1 128 V 7.314/ V °- 7811 / V 1-429/ V 7mJ y/ 

7c = 214.0 ccTdT (^L_y 3/r ( *_eIM Y 1/7 (j^v 11 ( _^_V /7 , (16) 

1 128 V 7.314/ V 0.7811 J V1.429/ V 7m Jy/ 

All these values are consistent with the results of previous studies (Liu et al. 2004, 2006). The 
electron-electron collision time, r c = 7 mgC 3 /87me 4 In A ~ 1164 72% 1 hrs (In A = 40 in this 



3 Here Z^Turb could be a complicated function of v s and v\- We, however, focus on the simplest case. 
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case), is much longer than any other time scale. Particle-particle collisions, therefore, may 
be ignored. The electron pressure (~ n^ c m e c 2 ) is about 5 times higher than the magnetic 
field pressure. The total particle pressure could be much higher. 

From Kirchhoff's law, the synchrotron absorption coefficient 

k u = E v j2^ c m e v 2 = (ire n/3\ / 3^B)[I(x M )/x M ] , 
we have the optical depth through the emission region 

t v = k u R = 2.482 C\ R\ 2 n 6 7 B^[I(x M ) /x M ] 

= 1.783 x 10- 7 Cr 5 C 2 - 16/7 D 8 8/7 (^L_Y 2 \ x mI(x m ) A 

1 2 8 V7.314/ |_ 0.7811 J ' V ; 

so the source is optically thin above x M = 1.766 x 1CT 3 , i.e., for v exceeding 34.52 GHz. 
Then the photon energy density 



12/7 n ~6/7 



U ~ £syn/4vrc R 2 = 127.9C^C~ LZ/ 'D 



X 



KtmV l^TSirJ Il429) \J^Tj) 6rgSCm (18) 

is more than 3 times lower than that of the magnetic field, justifying our neglect of the 
self-Comptonization cooling 4 . 

The self-Comptonization flux density is then given by (Blumenthal & Gould 1970) 
F x {v) ~ (2ire 7 n 2 BR 4 /3V3m 3 e c 6 D 2 )G(v/4v cl 2 ) = 2.121 n 2 B 1 R\ 2 DfG(v/Av cl 2 c )^y , (19) 
where 

roc pi 

G(z) = z dx dy exp(-x) I (zx^y' 1 ) (2 In y + 1 - 2y + y" 1 ) , (20) 
Jo Jo 

the low limit of the integral over x = r )/ r ) c has been extended from \/ r y c to 0, which simplifies 
the formula and introduces a less than 0.2% error, and z = iy/Ai/ c 'j 2 . Similarly one can define 
the spectral index 

a x = d \n(uF x )/d In v = 1 + d In G(z)/d In z . (21) 
The normalized power spectrum zG(z) and a x are depicted by thin lines in Figure 1 (left). 



4 The self-Comptonization effects may be incorporated into the model by replacing the B 2 with B 2 + 8nU 
in the electron cooling and source luminosity terms. 
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For the above mentioned NIR flare, the model predicts that in the Chandra and XMM 
bands (u x = z/j|l0 18 Hz, which corresponds to ~ 4.2 keV) 

z = 0.28 C^C-^D-^ 4 13/7 ( ZMlM) 2A (J*±.) " 13/7 ( -J}-) " (22) 

i 2 8 18 V7.3147 V 0.7811 y V 1.429/ V7mJ y y V ' ; 

ax(z) = 1.1, (23) 



F x {z) = 0.0033 C^C-^D-^ x 



/iJy, (24) 



.7.314/ V 0.7811 y V 1.429/ V7mJ y y 

which is about 5 times below the quiescent flux level. We therefore do not expect a detectable 
X-ray flux accompanying this flare. From equations (22) and (24), one obtains the correlation 
between the NIR and X-ray spectroscopic observations. For a given X-ray spectral index, 
which is a function of z, Equation (22) gives the dependence of the NIR flux density on 
xmi which determines a. For a given X-ray flux density F x , one obtains function 
of z by eliminating the F v terms in the two equations. This gives F v as a function of a 
with xm as a parametric variable. These two relations are plotted in Figure 1 (right) for 
C± = C 2 = D$ = 1.0. Here F x is in units of Sagittarius A*'s quiescent flux density of 
0.015 /iJy at 4.2 keV. 



5. Discussion and Conclusions 

First, we notice that for weak NIR flares with very hard spectra (to the bottom-right of 
the figure), there is little X-ray emission and that X-ray flares are expected to accompany 
NIR flares with soft spectra (to the left). These predictions may be readily tested via a 
statistical study of all the observed flare events (Hornstein et al. 2002) and explain why the 
occurrence rate of NIR flares is more than two times greater than that of the X-ray events. 
For NIR flares with very hard spectra, a x is also large (corresponding to a hard spectrum), 
we expect hard X-ray/7-ray emission, especially during bright NIR flares (to the top-right). 

Second, brighter X-ray flares (to the top-left) tend to have higher NIR fluxes and softer 
NIR spectra. a x is also small to the left. We expect little 7-ray emission accompanying 
weak NIR flares with very soft spectra (to the bottom-left). 

Third, for a given a, F x is roughly proportional to F v . The spectral index a x correlates 
with F v for low values of a, but this correlation almost vanishes for positive a's. For low 
values of a, we also expect a correlation between F x and a x . 

Finally, the results are very sensitive to the values of C\ and C 2 . The dashed lines in 
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the figure indicates X-ray flares with C\ = 0.5. These are quite different from the lines with 
C\ = 1.0. For the referenced NIR flare, the X-ray flux will be detectable with Chandra for 
C\ = 0.5. Thus C\ and C 2 may be determined accurately with SSO. 

Some of these predictions are shown explicitly in the left panel of Figure 2, where the 
solid line is the model spectra of the synchrotron and self-Comptonization components for the 
NIR flare studied above. Although there is little X-ray emission, we expect 7-ray emission 
up to a few MeV during the flare. The model is also applied to the brightest flare observed 
simultaneously in the NIR and X-ray bands by Eckart et al. during their multi-wavelength 
campaign in 2004. The NIR and X-ray fluxes correlated well during the flare and the peak 
flux densities F v = 7mJy and Fx = 0.223 /iJy. The spectral indexes of the flare have not 
been published yet. Because the X-ray flux is high, we expect soft NIR spectra. The dotted 
(with d = C 2 = 1.0) and dashed (with C x = 2.2, C 2 = 0.74) lines are for this flare. The 
right panel shows how the NIR and X-ray spectral indexes depend on C\ and C 2 . The cross 
and square correspond to the dotted and dashed lines in the left panel, respectively. These 
show that for given flux densities the spectral indexes are very sensitive to the values of C\ 
and C 2 - Because both C\ and C 2 are dimensionless quantities, we expect them to be on the 
order of 1. The model therefore predicts that the NIR emission has a soft spectrum with 
a < 0. 

For the referenced NIR flare the electron acceleration time 

is shorter than the dynamical time of r dy ~ 20 mins, which is consistent with the observed 
flare variation time, justifying our usage of the steady-state solution, which, however, may 
not be valid when r ac > Td y , as indicated by the shaded area in Figure 1 (right). Note that 
for the NIR flares with softer spectra we studied before (Liu et al. 2006), a < —2.0, the 
steady-state solution therefore may not be applicable. In the most general cases, both r ac 
and To are functions of t, which are determined by the large scale MHD processes. Time de- 
pendent solutions of equation (1) N(j,t) are needed. However, we don't expect the electron 
distribution to be radically different from Maxwellian. By introducing a time-dependent 
7 c (t), the results here may be readily generalized to incorporate the time evolution of the 
system and applied to the accretion flow in the quiescent state when advection can be im- 
portant. One may also have to take into account the Doppler and gravitational effects as 
the emitting plasma is likely moving on a Keplerian orbit near the black hole (Bromley et 
al. 2001; Broderick & Loeb 2005). These are beyond the scope of the current investigation. 

Since the discovery of NIR and X-ray flares in Sagittarius A*, several models have been 
proposed (Markoff et al. 2001; Liu & Melia 2002). Given the uncertainties in the electron 
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acceleration mechanism, earlier models prescribe either the electron distribution (Yuan et 
al. 2004) or the turbulence spectrum (Liu et al. 2004) with several free parameters re- 
sulting in all kinds of emission spectra. In a previous work (Liu et al. 2006), we showed 
that Sagittarius A*'s flare activity might be explained with thermal synchrotron and self- 
Comptonization emission and all the model parameters, namely B, n, R, and 7 C , could be 
readily determined by SSO. In this paper we have shown that the electrons are likely heated 
by turbulence, which produces a Maxwellian distribution and predicts two time- insensitive 
constraints on the plasma properties making the model testable with SSO. This is the first 
model that can be tested with current instrument capability. We have provided general 
expressions for the emission spectra and demonstrated how properties of the flaring plasma 
may be inferred from the spectroscopic observations or simultaneous flux measurements in 
both bands. The model may therefore provide powerful tools for studying physical processes 
near the event horizon of black holes. A time-dependent treatment of this problem will make 
it possible to compare MHD simulations of the black hole accretion with the observed flare 
activity in Sagittarius A*, which will eventually lead to a measurement of the black hole's 
spin after the instability is identified and light propagation effects are incorporated properly. 
(Note that for the referenced NIR flare the electron Lorentz factor 7 C ~ 200 already suggests 
that the gravitational energy release beyond the last stable orbit of a Schwarzschild black 
hole is not sufficient to power the flare.) The formulae for the thermal synchrotron and 
self-Comptonization emission are valid so long as optically thin synchrotron emission dom- 
inates the radiative output and may therefore be applicable to other astrophysical sources, 
particularly Low-Luminosity- AGNs. 
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Fig. 1. — Left: Thermal synchrotron (thick) and self-Comptonization (thin) power spectral 
index (dotted) and normalized energy flux density vF v (solid) as functions of their normalized 
frequencies. Right: Correlation between the NIR (1/14 = 1.429) and X-ray (vis = 1.0) spectral 
indexes and flux densities for C\ = C2 = D 8 = 1.0. The horizontal and vertical axises are 
for the NIR spectral index and flux density, respectively. The dotted and solid lines indicate 
constant X-ray spectral index and flux density (in units of the quiescent level of 0.015/iJy) 
with the corresponding value labeled on the lines. The dashed lines are similar to the solid 
lines but with C\ = 0.5. The dot indicates the model studied in the paper. The steady-state 
solution may not be valid in the shaded region. See text for details. 
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Fig. 2. — Left: Model fit to the NIR flare studied in the paper (solid line) and the peak 
flux densities of the brightest flare observed simultaneously in the NIR and X-ray bands by 
Eckart et al. (2006). For the latter, F u = 7 mJy and F x = 0.223 /iJy. C x = C 2 = 1.0 
for the dotted line, and C\ = 2.2, C 2 = 0.74 for the dashed line. Note that the effect of 
self-absorption is not included here. The Right panel shows the dependence of the model 
predicted NIR (a) and X-ray (ax) spectral indexes on G\ and C 2 for this flare. We expect 
a soft NIR spectrum with a < 0.0 for reasonable values of C\ and C 2 . 



